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ABSTRACT 

Preliminary  results  are  presented  from  a  numeri¬ 
cal  study  of  the  unsteady  hydrodynamics  of  an  AUV 
thruster.  A  vortex-lattice,  lifting-surface  model,  devel¬ 
oped  originally  for  ship  propellers  operating  under  nor¬ 
mal  design  conditions  at  constant  angular  velocity,  is 
modified  to  handle  unsteady  operating  conditions  of  an 
AUV  thruster  during  dynamic  positioning  and  maneu¬ 
vering.  The  present  numerical  code  applies  to  an  un¬ 
ducted  propeller  though  the  method  can  be  modified  to 
handle  ducted  propellers.  The  results  from  sample  runs 
of  the  thruster  undergoing  step  changes  and  sinusoidal 
oscillations  in  the  angular  velocity  of  the  propeller  are 
shown. 
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INTRODUCTION 

The  motivation  of  this  research  is  to  develop  an  accurate  h}'dro- 
d3Uiamic  model  of  an  AUV  thruster  during  unsteady  oi^crations. 
Such  a  model  could  be  used  to  develop  controllers  for  vehicles 
that  operate  in  the  ocean  and  are  subjected  to  large  time- varying 
forces.  For  example,  the  large  oscillating  wave  forces  that  dom¬ 
inate  near-shore  and  near-surface  ocean  environments  will  cause 
an  underwater  vehicle  to  move  back  and  forth  with  horizontal  and 
vertical  amplitudes  of  the  order  of  meters  and  with  time  scales 
of  the  order  of  seconds.  To  maintain  a  prescribed  course  under 
closed-loop  control,  the  thrusters  will  have  to  continually  switcli 
from  forw^ard  to  reverse  thrust  and  back  again  to  forward  thrust 
in  such  a  way  as  to  negate  the  effects  of  wave  and  current  forces  as 
well  as  the  inertial  forces  of  the  vehicle.  The  correct  timing  of  this 
control  action  is  possible  only  if  there  is  an  accurate  model  of  the 
thruster  dynamics  and  hydrodynamics  built  into  the  controller. 

When  a  d3uiamic  positioning  controller  commands  a  certain 


amount  of  thrust  using  present-da}’  methods,  the  desired  tlirust 
does  not  occur  instantaneously  becau.se  of  time  delay’s  involving 
the  fluid  inertia  and  velocity  of  the  fluid  flowing  over  the  pro¬ 
peller  blades.  Consider  a  propeller  operating  with  a  steady-state 
angular  velocity  0  and  forward  thrust  T.  If  the  angular  velocity 
is  suddenh'  reversed  to  —H,  the  fluid  surrounding  the  propeller 
will  still  be  flowing  in  the  original  direction.  With  time,  the  fluid 
reverses  course  and  evcnttially  reaches  a  steady-state  flow  velocity 
equal  to  the  original  value  but  in  the  opposite  direction  (a.ssuining 
that  the  thruster  is  S3’minetric  for  forward  and  reverse  operations). 
The  thrust  also  reverses  direction  and  becomes  equal  to  —T  at 
steady  state.  This  process  of  reversing  thrust  from  the  original 
command  to  the  actual  application  of  the  desired  value  can  take 
a  number  of  seconds  on  an  AUV,  enough  of  a  delay  to  cause  a 
characteristic  limit  cycle  which  reduces  tracking  performance  in 
dynamically  positioned  vehicles. 

\oerger  et  al.  (1990)  were  the  first  to  recognize  that  thruster 
d3mamics  were  responsible  for  limit  c3^cles  in  the  context  of  d}’- 
iiamic  positioning  of  underwater  vehicles.  They  proposed  a  first 
order  d3mamic  model  to  describe  this  phenomenon.  The  model 
includes  a  term  in  the  momentum  equation  that  is  linearly  re¬ 
lated  to  the  angular  acceleration  of  the  motor  shaft  as  well  as  the 
square  of  the  shaft  angular  velocity.  The  inertia  coefficient,  which 
is  found  from  experiments,  implicitly  accounts  for  the  added  mass 
of  the  fluid.  However,  this  model  ignores  the  effect  that  fluid  flow 
has  on  the  instantaneous  angle-of-attack. 

Healey  et  al.  (1994)  derived  a  second  order  model  that  uses 
the  angular  velocity  of  the  propeller  shaft  and  the  axial  fluid  ve¬ 
locity  as  state  variables.  This  formulation  permits  the  angle-of- 
attack  of  the  propeller  blade  to  vary  with  time.  Simulations  asing 
the  model  showed  that  the  prediction  of  both  steady-state  and 
transient  responses  was  improved.  Whitcomb  and  Yoerger  (1995) 
simplified  this  model  by  writing  it  as  a  one  dimensional  system 
with  the  axial  fluid  velocity  as  the  state  variable  and  the  angular 
velocit3"  of  the  propeller  as  the  control  input. 

The  accurate  real-time  measurement  of  axial  flow  velocity 
could  be  used  with  the  model  of  Whitcomb  and  Yoerger  (1995) 
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Figure  1:  Schematic  of  propeller  and  wake  with  distributed 
sources  and  vortices  (from  Kerwin  and  Lee  1978). 


Figure  2:  Numerical  simulation  of  development  of  vortex  wake 
sheet  behind  a  single  bladed  propeller  operating  with  variable 
angular  velocity.  The  wake  trails  from  the  i)ropeller  at  the  left 
to  downstream  on  the  right.  The  initial  calculations  (represented 
by  the  farthest  downstream  sections  of  the  helical  vortex  sheet) 
correspond  to  500  RPM.  The  angular  velocity  was  then  reduced 
to  100  RPM  causing  the  coils  of  the  wake  to  spread. 


to  design  a  thrust  controller.  However,  it  is  not  practical  to  mea¬ 
sure  the  fluid  velocity  near  the  propeller  blade  of  an  AUV  oper¬ 
ating  in  the  ocean.  This  is  where  the  numerical  hydrodynamic 
model  for  the  thruster  could  prove  useful.  It  would  be  used  to 
construct  an  observer  that  estimates  the  fluid  velocity  from  the 
past  history  of  the  propeller’s  angular  velocity  and  motor  torque 
(which  is  proportional  to  the  electric  current)  both  of  which  are 
easily  measured  in  real  time.  The  accuracy  of  the  observer,  then, 
depends  on  the  accuracy  of  the  hydrodynamic  model. 

DESCRIPTION  OF  NUMERICAL  MODEL 

The  numerical  model  is  based  on  the  vortex-lattice,  lifting-surface 
theories  developed  at  Massachusetts  Institute  of  Technology  and 
reported  in  Kerwin  and  Lee  (1978)  and  Keenan  (1989).  In  this 
formulation,  the  propeller  geometry  and  flow  field  are  represented 
by  a  set  of  sources  and  vortices  distributed  on  the  mean  camber 
surface  of  each  blade  and  a  distribution  of  vortices  shed  into  the 
wake  (Fig.  1).  The  strengths  of  the  sources  and  vortices  as  well 
as  the  location  of  the  shed  vortices  vary  with  time  and  can  be  de¬ 
termined  from  the  “no-flow”  boundary  condition  of  the  propeller 
(i.e.  the  normal  flow  velocity  is  zero  on  the  propeller  blades)  and 
the  “no-jump”  pressure  condition  on  the  wake  vortex  sheet  (i.e. 
the  fluid  pressure  is  continuous  across  the  vortex  sheet). 

Keenan  (1989)  developed  vortex-lattice  methods  to  simulate 
the  dynamics  of  a  thruster  subjected  to  a  spatially  varying  and 
time  invariant  in-flow  current.  In  this  scheme,  vorticity  is  slied 
at  each  time  step  from  the  tip  and  trailing  edge 'of  the  propeller 
into  the  wake.  The  strength  of  the  shed  vorticity  will  change  due 
to  the  varying  position  of  the  blade  and  in- flow  velocity.  Since, 
the  pattern  repeats  itself  every  revolution  of  the  blade,  an  steady- 
state  oscillatory  pattern  develops  in  the  wake. 

We  have  recently  extended  the  method  so  that  the  in- flow 
current  and  the  angular  velocity  of  the  propeller  can  vary  with 
time.  Fig.  2  shows  a  sample  calcidation  of  the  development  of 
the  wake  due  to  a  step  change  in  the  angular  velocity  from  500 


to  100  revolutions  per  minute  (RPM)  and  vehicle  velocity  of  2.5 
m/s. 

Our  time-variant  model  is  presently  limited  to  an  unducted 
propeller  and  none-reversing  flow.  Ducted  propellers  have  been 
simulated  in  spatially  varying  and  time  invariant  flows  using  vortex- 
lattice  methods  (Hughes  1993),  and  it  is  a  straight  forward  ex¬ 
tension  to  temporally  varying  the  in-flow  currents  and  propeller 
angular  velocity.  Reverse  flow  is  more  difflcult  to  imi^lcmeiit  be¬ 
cause  of  the  need  to  track  the  vortex  sheet  as  it  is  pulh'd  back 
through  the  propeller.  Our  future  work  will  involve  developing 
a  model  that  permits  each  vortex  segment  to  have  compact  suj)- 
port  and  convect  independently  (Anderson  and  Grecngard  1985). 
Though  this  scheme  would  fail  to  simulate  the  turl^ulence  that 
would  develop  during  the  “ingestion”  of  the  vortex  sheet,  it  will 
help  us  understand  the  large-scale  dynamics  of  the  jiroccss. 

Another  limitation  of  the  numerical  code,  in  its  present  form, 
is  its  ability  to  simulate  operating  conditions  that  are  away  from 
the  design  conditions  of  the  thruster  as  defined  by  the  advance 
coefficient  of  a  propeller  which  is  defined  as: 


.where  U  is  the  ambient  fluid  velocity  due  to  vehicle  motion  and 
current,  n  is  the  angular  velocity  of  the  propeller,  and  D  is  the 
diameter  of  the  propeller.  Typical  propellers  are  designed  to  op¬ 
erate  at  a  specific  value  of  J.  AUV  thrusters  on  the  other  hand 
operate  over  a  range  of  J-valucs.  Because  of  this,  large  aiiglcs  of 
attack  arc  possible  as  well  as  separation  about  the  leading  edge  of 
the  propeller  blades.  Further  modifications  of  the  numerical  code 
will  include  this  effect  through  the  modification  of  the  leading- 
edge  suction  force  (Keenan,  1989). 

SAMPLE  CALCULATIONS 

In  this  section,  we  present  results  from  three  separate  simulations. 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
0MB  No.  0704^788 


ll'#gORTDATE  3.  REPORT  TYPE  AND  OATES  COVERED - - 

4.  TITLE  ANO'SUBTITLE  - ^  - - — I _ TECHNICAL 

NUMPRTrar  a.  i-unding  NUMBERS  ^ 

_ '  ONR  N00014-96-15014 

6.  AUTHOR(S)  - - - — _ _ _ _ _ 

J.  KNOWLES,  M.  GROSENBAUGH  AND  D.  KEENAN 


7.  PERFORMING  ORGANIZATION  NAME{S)  AND  AOORESS(ES) - - - 

WOODS  HOLE  OCEANOGRAPHIC  INSTITUTION 
WOODS  HOLE,  MA  02543 

».  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADORESS(ES) - 


ORGANIZATION 
REPORT  NUMBER 


WHOI  CONTR. 


MONITORING 

AGENCY  REPORT  NUMBER 


OFFICE  OF  NAVAL  RESEARCH 
ENVIRONMENTAL  SCIENCES 
DIRECTORATE 

ARLINGTON,  VA  22217-5660 

11.  SUPPLEMENTARY  NOTES - - - - - -  I 

In  atmg  this  report  in  a  bibliography,  the  reference  given  should  be: 

ISOp'l^iLf  “■>  conference 

12a.  DISTRiHimnM/AtfiiM 

lib.  DISTRIBUTION  CODE 

approved  for  PUBLIC  RELEASE: 

distribution  unlimited 


13.  ABSTRACT  (Maximum  200  wortisf 


Preliminary  results  are  presented  from  a  numeri- 
ca  s  udy  of  the  unsteady  hydrodynamics  of  an  AUV 
thriKter.  A  vortex-lattice,  lifting-surface  model,  devel¬ 
oped  originally  for  ship  propellers  operating  under  nor- 
mal  design  conditions  at  constant  angular  velocity,  is 
mo^fied  to  handle  unsteady  operating  conditions  of  an 
AUV  thruster  during  dynamic  positioning  and  maneu¬ 
vering  The  present  numerical  code  applies  to  an  un- 
ducted  propeller  though  the  method  can  be  modifled  to 
handle  ducted  propellers.  The  results  from  sample  runs 
of  the  thruster  undergoing  step  changes  and  sinusoidal 
shown**””*  angular  velocity  of  the  propeller  are 


14.  SUBJEa  TERMS 


1)  AUV  thruster 

2)  unsteady  propeller 

3)  hydrodynamics 


15.  NUMBER  OF  PAGES 

4 

16.  PRICE  CODE 


°^UNCLAS<;iFnrTr*'  20.  UMITATIONOFABSTRAa 

UNCLASSIFIED  UNCLASSIFIED  UNCLASsfe® 


NSN  7540-01-280-5500 


Number  of  Vorlex  Panels  Retained  In  Wake 


Time  (sec) 


Figure  3:  Convergence  of  steacl3'-state  solution  as  a  function  of 
number  of  vortex  panels  in  the  wake.  The  angular  velocity  of  the 
propeller  is  700  RPM  and  the  flow  velocity  is  1.5  in/s. 

The  first  is  a  propeller  operating  under  steady  conditions.  The 
second  simulation  is  a  propeller  undergoing  sinusoidal  variations 
of  its  angular  velocity  about  a  steady-state  operating  condition. 
The  final  case  is  a  step  test  where  the  angular  velocity  and  flow 
velocity  is  increased  instantaneously  from  zero  to  a  finite  value. 

The  propeller  that  we  model  is  a  three-bladed  propeller  with 
a  diameter  of  24  cm  manufactured  by  Vetus  Corporation  (model 
number  BP125).  This  is  the  same  propeller  that  is  used  on  the 
JASON  ROV,  and  it  was  used  in  the  experiments  of  Whitcomb 
and  Yoerger  (1996). 

Steady-State 

The  stead\'“State  case  is  meant  to  verify  the  stability  of  the  code 
and  to  determine  at  what  point  the  wake  should  be  truncated 
downstream.  The  simulation  begins  by  initially  specifying  an 
ambient  flow  velocity  [/,  the  propeller’s  angular  velocity  Q,  the 
number  of  vortex  panels  on  the  blades,  and  the  number  of  pan¬ 
els  that  make  up  the  wake.  Next,  we  calculate  the  steady-state 
solution  using  a  steady-state  code  by  adjusting  the  strength  of 
the  vortex  panels  on  the  blade  and  in  the  wake  in  order  to  satisfy 
the  “no-flow”  boundary  condition  on  the  blade  and  the  “no-flow” 
condition  of  the  wake  (Kerwin  and  Lee,  1978).  Finali}^  we  take 
this  solution  as  the  starting  point  of  a  time-domain  simulation. 
We  move  the  propeller  a  fraction  of  a  revolution  (in  this  case  12'’), 
shed  the  row  of  vortex  panels  off  the  trailing  edge  into  the  ncar 
wake,  convect  them  downstream,  and  remove  the  row  of  vortex 
panels  from  the  far-end  of  the  wake.  If  our  time-domain  code 
is  stable,  then  the  steady  state  solution  will  persist  for  all  time, 
which  it  does. 

Fig.  1  presents  the  results  for  the  case  of  f/  =  1.5  m/s  and 
n  =  700  revolutions  per  minute  (RPM).  The  length  of  the  simu¬ 
lation  is  1.5  s.  The  number  of  panels  used  in  the  wake  is  varied  in 
order  to  determine  when  the  code  converges.  Since  truncating  the 
wake  is  an  approximation,  the  length  of  the  wake  will  effect  the 
accuracy  of  the  calculations.  Here  convergence  is  achieved  with 
2400  panels  based  on  the  values  of  the  thrust  and  torque. 


Figure  4:  Thrust  and  torque  of  a  propeller  undergoing  simi.soidal 
variations  of  the  propeller’s  angular  velocity  about  a  mean  oper¬ 
ating  condition.  Tlie  mean  angular  velocity  is  700  RPiM  and  the 
amplitude  of  the  variation  is  ±200  RPM.  The  flow  velocity  is  1.5 
m/s. 

Sinusoidal  Variation  of  Propeller  RPM 

For  this  case,  the  initial  conditions  are  the  same  as  for  the  pre¬ 
vious  steady-state  simulation  with  a  wake  length  of  2400  panels. 
The  simulation  proceeds  as  before  for  one  proj^eller  cycle  after 
which  the  value  of  Q  is  varied  sinusoidally  from  500  RPM  to  900 
RPM  with  a  period  of  0.5  seconds  which  is  one-half  of  a  propeller 
revolution  at  the  nominal  angular  velocity  of  700  RPM.  The  re¬ 
sults  of  the  simulation  show  a  time  lag  between  the  peaks  of  the 
angular  velocity  and  the  peaks  of  the  thrust  and  torque  curves 
(Fig.  4).  The  mean  values  of  the  thrust  and  torque  match  the 
steady-state  values. 

Step  Test 

The  step  test  involves  an  instantaneous  jump  from  =  0  and 
C/  =  0  at  time  i  =  0  to  Q  =  700  RPM  and  U  =  m/s.  The 
initial  shape  of  the  wake  is  taken  as  the  steady-state  solution  with 
2400  vortex  panels,  but  the  initial  strengths  of  the  vortex  panels 
are  set  equal  to  zero.  During  the  first  time  step,  the  first  row  of 
vortex  panels  that  is  shed  into  the  wake  is  strong.  This  starting 
vortex  influences  the  solution  while  it  is  in  the  near  wake,  lo.scs 
influence  as  it  convccts  downstream,  and  has  no  influence  once  it 
drops  off"  the  end  of  the  truncated  w'ake.  This  is  seen  in  the  thrust 
and  torque  results  (Fig.  4). 

The  peaks  in  the  thrust  and  torque  curves  are  due  to  the 
influence  of  the  starting  vortex  which  is  shed  at  time  t  =  0.  The 
first  peak  corresponds  to  roughly  one-third  of  a  propeller  revo¬ 
lution  and  corresponds  to  the  three  propeller  blades  j^assing  the 
starting  vortices  shed  by  the  blade  that  is  rotating  ahead  of  it. 
The  next  peak  at  roughly  two- thirds  of  a  revolution  corrcsi)ouds 
to  the  trailing  blades  passing  the  starting  vortices  which  have 
convected  further  downstream.  A  smaller  peak  at  one  proi)eIler 
revolution  corresponds  to  the  blades  pa.ssing  their  own  starting 
vortices.  The  thrust  and  torque  values  asymtote  to  the  steady- 
state  solution. 


Figure  5:  Thrust  and  torque  of  a  propeller  undergoing  an  instan¬ 
taneous  acceleration  from  0  RPM  to  /^OO  RPM.  The  flo^v  velocitj'^ 
is  assumed  to  go  instantly  from  0  m/s  to  L5m/s. 

SUMMARY  AND  CONCLUSIONS 

We  have  begun  a  research  program  to  simulate  the  hydrodynam¬ 
ics  of  an  unsteady  AUV  thruster  with  the  unsteadiness  arising 
from  temporal  changes  in  the  inflow  velocity  and  the  propeller’s 
angular  velocity.  The  present  code  applies  to  an  unducted  pro¬ 
peller  operating  in  a  regime  in  which  the  vortex  sheet  is  always 
convected  downstream  in  one  direction.  Future  improvements 
will  include  the  ability  to  simulate  a  ducted  propeller  and  fully- 
reversed  flow  where  the  wake  is  drawn  back  past  the  propeller 
blades. 

The  motivation  for  developing  the  simulation  is  to  use  it  to 
design  closed-loop  thrust  controllers  for  AUVs.  The  results  of 


simulations  show  the  significant  effect  of  fluid  dynamics  on  the 
thrust  and  torque.  Because  torque  can  be  measured  through  its 
relation  with  motor  current,  our  results  indicate  that  torque  could 
be  used  to  construct  an  accurate  observer  for  the  flow  vc^Iocity  at 
the  propeller  blades  which  would  account  for  the  hydrodynamic 
effects  on  the  thrust. 
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